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ABSTRACT 

Based on recent findings of a form ation mechanism of substructure in tidal tails by 
iKupper. Macleod fc Hegglel ( 2008) we investigate a more comprehensive set of TV-body 
models of star clusters on orbits about a Milky- Way-like potential. We find that the 
predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster 
follows as long as the cluster lives long enough for the overdensities to build up. 

The distance of the overdensities along the tidal tail from the cluster centre de- 
pends for circular orbits only on the mass of the cluster and the strength of the tidal 
field, and therefore decreases monotonically with time, while for eccentric orbits the 
orbital motion influences the distance, causing a periodic compression and stretching 
of the tails and making the distance oscillate with time. We provide an approximation 
for estimating the distance of the overdensities in this case. 

We describe an additional type of overdensity which arises in extended tidal tails 
of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar 
stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or 
a disk shock is not the direct origin of an overdensity within a tidal tail. Escape due to 
such tidal perturbations does not take place immediately after the perturbation but 
is rather delayed and spread over the orbit of the cluster. All observable overdensities 
are therefore of the mentioned two types. In particular, we note that substructured 
tidal tails do not imply the existence of dark-matter sub-structures in the haloes of 
galaxies. 

Key words: galaxies: kinematics and dynamics - galaxies: star clusters - methods: 
analytical - methods: A^-body simulations - galaxies: haloes - cosmology: dark matter 



1 INTRODUCTION 

Since an increasing number of tidal st ructures are being 
disco v ered around the Milky Way (e.g . Odcnkirchcn et ahl 
200ll: iBelokuroy et all l200rj: iGrillmair fe Dionatosl l200rj; 
Grillmaidl2009l: lJuric et al.ll2008l ; iKeller. Da Costa fe Priori 
20091 : iKoch et al.ll2009D ." there is a profound interest for a 
deeper physical understanding of tidal tails. Tidal tails offer 
direct constraints on the ongoing dissolution process of star 
clusters as well as on the last few Gyr of their evolution. 
Furthermore, their alignment, shape and structure give a 
unique opportunity to reconstruct the cluster's orbit in un- 
precedented detail and thus probe our understanding of the 
Milky- Way potential. 

From the simple assumption that a star cluster has a 
nearly constant mass-loss rate on its orbit about the galactic 
centre, its tidal tails are often expected to have a more or less 
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homogenous and smooth structure. Substructure in the tidal 
tails is usually only expected in time-variable tidal fields 
when the mass-loss rate varies with time. 

However, in a recent paper (Kiipper, Macleod & Heg- 
gie 2008; hereafter KMH) we analytically predicted and nu- 
merically proved the formation of substructure in tidal tails 
even within a constant tidal field. This effect is due to an 
epicyclic motion of stars evaporating from the cluster in- 
stead of the more common approximation of linear motion 
along the tails (for details see KMH). Like a standing wave 
this motion leads to statistical overdensities if it is performed 
by a stream of stars. 

We were able to show that these epicyclic overdensities 
arise at a distance of many tidal radii along the orbit from 
the cluster centre and also gave a simple relation between 
this distance and the tidal radius of the corresponding clus- 
ter. The numerical verification was done for a fairly simple 
iV-body model, using a point-mass galaxy and a star clus- 
ter without a stellar-mass spectrum, to avoid all disturbing 
influences. 

More recently, this effect has also been studied by 
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lJust et all (|2009h in a more general context, showing that, 
even in a Milky- Way-like potential and with a cluster con- 
sisting of stars drawn from a mass function, these over- and 
underdensities arise and can give valuable information on 
the current state of the cluster. 

In KMH we found that the prerequisite for the forma- 
tion of epicyclic overdensities is that a majority of escaping 
stars leaves the cluster with velocities slightly above the es- 
cape velocity through one of the two Lagrange points. If the 
scatter in escape conditions is too large, the statistical over- 
densities are expected to vanish. For constant tidal fields 
this condition has been shown to be well fulfilled, but for 
time- variable tidal fields this has not been investigated yet. 

For star clusters on eccentric orbits or on orbits in- 
volving periodic disk shocks the tidal field varies with time 
and the internal evolution of the cluster is perturbed. Thus, 
the scatter in escape conditions may be increased such that 
epicyclic overdensities may not be able to arise. How large 
the influence of the perturbation is depends on the orbital 
parameters. Therefore, we here investigate a comprehensive 
set of TV-body simulations of globular clusters on different 
orbits about a Milky- Way-like potential, and study the for- 
mation of substructure within the tidal tails (Sec. [3}. Within 
the models we also look for overdensities in the tails which 
can be assigned to a tidal variation and an associated change 
in the mass-loss rate. But first a few theoretical considera- 
tions have to be set out to explain the various aspects of es- 
cape and, in this context, explain our methodology (Sec. [2}. 



2 THEORY OF ESCAPE 

Understanding the formation of tidal tails, and especially 
of substructure within those tails, is relatively easy in the 
ideal case of a star cluster on a circular orbit about a spher- 
ically symmetric galaxy as studied in KMH. Escape from 
such a cluster can be theoretically described in the frame- 
work of epicyclic theory in which the equations of motion 
of escaping stars are fairly simple. KMH showed that stars, 
when they evaporate with small velocities from the cluster 
through one of the Lagrange points, move on oscillatory or- 
bits along the tidal tails. These oscillations lead to periodic 
stellar over- and underdensities when they are performed by 
a constant stream of escapers. But how does this change for 
more complex clusters in time- variable tidal fields? 

Since the formation of substructure in the tails, as de- 
scribed in KMH, is due to stars evaporating from the cluster, 
and in a constant tidal field evaporation is mass-loss driven 
by t wo-body relaxation l)Kiipper. Kroupa fc Baumgardt! 
120081 ). we argue that the more stars escape with evapora- 
tive escape conditions from the cluster, i.e. escape with low 
speed through one of the Lagrange points, the more pro- 
nounced are the overdensities and the easier it is to observe 
overdensities beyond the first one. 

In time-variable tidal fields this does not necessar- 
ily mean that a) escape has to be due to two-body re- 
laxation, since pericentre passages during eccentric orbits 
and disk shocks also create e scapers, as they enhance 
the mass-loss rate in gen eral l)Vesperini fc He'ggid Il997l : 
iBaumgardt fc M akino 20031) or that b) escape has to happen 
from the Lagrange points, as the distance of the Lagrange 



points from the cluster centre may vary more quickly on 
eccentric orbits than the cluster profile does. 

On the other hand, time variations in the tidal field, and 
associated variations in the mass-loss rate, are also the most 
frequently suggeste d formation m echanism for substructure 
in tidal tails i Dehnen et al]l2004h . However, this mechanism 
has not been convincingly proven in numerical experiments 
yet. Therefore, we need to study the mass- loss rate and the 
escape conditions of the investigated models in detail to un- 
derstand the origin of detected substructure. 

2.1 Mass-loss rate 

2.1.1 Constant tidal fields 

The time scale of evaporation in constant tidal fields is the 
relaxation time, t r h, as the mass loss of a tidally limited star 
cluster can be approximated by 

dM 



3/4 



(1) 



l|Baumgardtll200ll ). T he median two -body relaxation time : 
furthermore given by (|Spitzej|l987h 



1/2 R 3/2 



t rh = 0.138- 



'GV2 m i/2mA' ^ 
where N is the number of stars, Rh is the half-mass radius, 
G the gravitational constant, m the mean mass (M/N) and 
In A the Coulomb logarit hm where A is of order 0.11N in 
the case of equal masses (jGiersz fc Hcggic 1994a). 

For clusters on circular orbits the tidal field is static 
and the cluster will dissolve as a result of two-body relax- 
ation and the interactions of stars with the tidal field, where 
the mass-loss rate is determined appr oximately by the ra- 
tio of half-mass radius to tidal radius ()Gieles fc Baumgardt! 
2008). For such time-independent tidal fields the tidal radius 
scales with M 1 ^ 3 , hence decreases rather slowly with ongo- 
ing mass loss. Since the change in the half-mass radius also 
happens rather slowly, the mass-loss rate is approximately 
constant during one revolution about the galactic centre. 
Hence, substructure in the tidal tails can only arise through 
the epicyclic motion of escaping stars. 



2.1.2 Time-variable tidal fields 

If the cluster moves on an eccentric orbit or if the cluster 
periodically crosses the galactic disk then additional energy 
is put into the cluster on an orbital time scale and the mass- 
loss rate is expected to vary during one revolution. 

As the galactocentric distance of the cluster periodically 
changes on an eccentric orbit, the tidal radius is not mono- 
tonically decreasing any more but varies during one period. 
Furthermore, because the structure of the cluster may not 
be able to change rapidly enough, a part of the cluster may 
(temporarily) lie beyond the tidal radius when the cluster is 
at pericentre, where the tidal radius reaches its minimum. 
This induces a strong pull on those outlying stars, increas- 
ing on average the stellar energies. When the tidal radius in- 
creases again, a fraction of these once outlying stars may be 
recaptured but now with increased energy. The rest will be 
lost from the cluster. This whole process is often loosely re- 
ferred to as bulge shocki ng (for a detailed discussion of tidal 
shocking see for example iGnedin. Lee fc Ostrikerlll999l ). 
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During a disk shock the cluster is temporarily com- 
pressed on a time scale smaller than the internal dynamical 
time scale. As a consequence of this, all stars get a kick, 
accelerating some stars (mostly the stars in the outer parts) 
to velocities above the escape velocity, and these are then 
likely to leave the cluster. This behaviour can be described, 
for example, in the formalism of th e first- and second-order 
theory of disk shocks as reviewed in lGnedin. Lee fc Ostrikerl 
(1999). Here, the mean energy changes, (AE) and (AE 2 ), 
of the stars within the cluster can be expressed as 



(AE) 
{AE 2 ) 



2 2 

<x g m r , 

2 2 2 

« g m v r , 



(3) 
(4) 



where g m is the maximum vertical gravitational acceleration 
of the disk while r and v are the positions and velocities of 
the corresponding stars. The dependence on r 2 makes clear 
that, on the one hand a cluster is more affected by a shock 
when it is more extended, and on the other hand that the 
outer stars are most likely to be pushed to energies above the 
escape energy. Moreover, the energy shifts from equations [3] 
and [4] become less influential with increasing galactocentric 
radius of the disk sh ock since the vertical grav itational ac- 
celeration decreases. IVesperini fc Heggiel (|l997h found that 
disk shocks at solar radius and beyond are almost negligible 
for realistic globular cluster configurations. 

2.1.3 Delayed escape 

In both types of time-dependent tidal-field, pericentre pas- 
sages and disk shocks, the stars most affected are at radii 
comparable to the (apogalactic) tidal radius. At these radii 
the orbital time scale of stars within the cluster is on the or- 
der of the orbital time scale of t he cluster within the g alaxy . 

For constant tidal fields iFukushige fc Heggiel (|2000| ) 
have shown that a star which is pushed to an energy just 
slightly above the escape energy cannot escape from the clus- 
ter easily but has to pass near one of the two Lagrange 
points, as the potential barrier is lowest there. This sig- 
nificantly increase s the escape time of stars from a cluster 
|Baumgardtll200ll ). 

What happens in the case of time-dependent tidal fields 
has not been investigated yet. Nevertheless, for somewhat 
non-circular orbits it is likely that, for many affected stars, 
escape will be delayed until well after the pericentre passage 
or the disk shock. This delay in escape of a specific star may 
even last several revolutions, since the star may need several 
orbits within the cluster to escape, as in the circular case. 
Hence the time variation of the tidal field need not lead to 
de tectable overdensities within the tidal tails as suggested 
by iDehnen et alj (120041 ) if the fraction of delayed escapers 
is large and escape events are spread in time. However, the 
formation of epicyclic overdensites will still be possible, even 
if a majority of stars leaves the cluster as a result of tidal 
variations, as long as the escape conditions of the escaping 
stars are evaporation-like, i.e. escape with low velocities. 



2.2 Escape conditions 

2.2.1 Constant tidal fields 



Gicrsz fc Heggiel (ll994bh.lBaumgardt. Hut fc Heggiel (2002) 
and iKiipper. Kroupa fc Baumgardt ( 20081 ) were able to 



show that for clusters (in constant tidal fields) the escape ve- 
locities of the majority of escaping stars follow a log-normal 
distribution when normalized to the specific velocity disper- 
sion of the cluster. 

This behaviour can be understood in terms of two-body 
relaxation: the velocity distribution within the cluster is al- 
most Maxwellian and stars in the high-velocity tail of this 
distribution are deemed to leave the cluster as their veloc- 
ities exceed the escape velocity of the cluster. The latter is 
furthermore linked to the velocity dispersion for clusters in 
virial equilibrium. Hence, the ratio of the escape velocity of 
an escaping star to the velocity dispersion of the cluster at 
that time is expected to be constant. And since the velocity 
dispersion changes rather slowly with time, stars in the tidal 
tails have quite similar escape conditions (which we call here 
evaporative escape conditions) leading to thin and dynami- 
cally cold tidal tails with observable epicyclic overdensities. 



2. 2. 2 Time-variable tidal fields 

As mentioned above, tidal variations affect the velocities 
of cluster stars. Hence, for strong tidal variations the ve- 
locity distribution does not necessarily have to be nearly 
Maxwellian any more. Furthermore, the tidal radius and 
thus the escape velocity of the cluster changes during one 
revolution about the galactic centre. In addition, the veloc- 
ity dispersion periodically varies with time as the cluster 
gets compressed or oscillates. Finally, the cluster can even 
be temporarily out of virial equilibrium. 

From these considerations we may expect that the es- 
cape conditions have a larger spread around a mean value 
and may even deviate from log-normal and become asym- 
metric. In the extreme case of a cluster which gets com- 
pletely disrupted through a single shock the scatter in escape 
conditions will be maximal. Thus, broad and dynamically 
hot tidal tails without epicyclic overdensities are expected 
from this crude picture. The intermediate regime, from per- 
fectly evaporative escape conditions to maximally perturbed 
conditions, will be the subject of our parameter-space study. 



2.3 Number- and orbital- velocity distribution 

We are going to analyse the tidal tails of our models in 
terms of number- and orbital-velocity distributions, in or- 
der to try to find the origin of any observed overdensi- 
ties and relate them to the tidal radius of the cluster. In 
other words, motivated by the results of KMH, we investi- 
gate the number distribution of stars along the tidal tails 
and compare it to the average orbital velocities of stars at 
that position within the tails. Epicyclic motion of a con- 
stant stream of stars leaving the cluster creates over- and 
underdensities at periodic distance intervals. For constant 
tidal fields the distance of the density enhancements only 
changes when the conditions of escape are modified, e.g. 
when the tidal radius changes. Furthermore, the oscillatory 
(epicyclic) movement also creates a signal in the orbital ve- 
locity of the stars along the tails, as was f i rst di scovered by 
ICapuzzo Dolcetta. Di Matteo. fc Miocchll ([2005) in numeri- 
cal simulations: places where the stars move fastest (slowest) 
should correspond to the positions of underdensity (overden- 
sity) (cf. KMH). 



4 A.H.W. Kiipper, P. Kroupa, H. Baumgardt and B.C. Reggie 



/V(Rj) 




Figure 1. Sketch of the analysis setup. The grey ellipse represents 
the cluster while the small grey circle shows a star in the tidal tails 
with index i on an epicyclic motion within the tails. The solid line 
with black knots gives the orbit of the cluster, where the knots 
show equidistant points along this orbit. Rj is the galactocentric 
radius of the knot with index j, and V (Rj) the orbital velocity 
of the cluster at this position respectively. As star i is closest to 
the knot at Rj it is given the index ij and counted in bin j. Its 
velocity is Vij. is the tidal radius, and yc gives the length of 
the epicyclic loops which are indicated by the long-dashed lines. 



To analyse the clusters in a consistent way, all stars are 
binned along the orbit in bins of size 25 pc. For this purpose, 
we calculate points, Rj, along the orbit of the cluster at 
separations of 25 pc and assign the tail stars to the point 
which they are closest to. The resulting distribution gives 
the number density of stars as a function of distance from 
the cluster centre along the orbit. In addition, the average 
orbital velocity with respect to the cluster within the jth 
bin, Vj, is measured as 



\ i=0 



(5) 



where My is the speed of the ith star in the jth bin, Nj 
the number of stars in the jth bin and V (Rj) is the speed 
the cluster will have or had at the position Rj of the orbit 
(Fig. Furthermore, the velocity dispersion in each bin is 
calculated with 



(6) 



where Vj is the mean velocity in the jth bin and w| the mean 
squared velocity respectively. 



2.4 Tidal radius and distance yc 

In KMH we theoretically predicted and numerically verified 
a relation between the positions of the overdensities, yc, and 
the tidal radius of the cluster, xl- For a cluster in a point- 
mass galactic potential the position of the first overdensity 



yc = ±12kxl, 



but for the general case the relation gets more complicated. 
The distance of the first overdensity is then given by 



yc 



AM 



4ft" 



(8) 



where ft is the angular velocity of the cluster on its orbit 
about the galactic centre and k is the so-called epicyclic 
frequency (see KMH). The latter comes from epicycle the- 
ory, where the motion of a star is approximated by the cir- 
cular motion of a guiding centre about the galactic centre 
and an oscillation of the star about this guiding centre (cf. 
iBinnev fc Tremaind[l987l ). 

In this framework the epicyclic frequency for a cluster 
moving in a planar orbit is given by 

* 2 = ^ +3n ' (9) 
Here, $ is the galactic potential and R is the galactocentric 
radius of the cluster. The epicyclic frequency gives a quan- 
tification of the steepness of the galactic potential in the ra- 
dial direction and the strength of the centrifugal force, and 
therefore naturally appears in the equation for the tidal ra- 
dius (as the distance of the Lagrange points from the cluster 
centre is usually called): 



GM 



4ft 2 



(10) 



where G is the gravitational constant. The epicyclic fre- 
quency obviously has a large influence on the value of xl- 
While the relation k = ft holds for point-mass galaxies, 
the value for a Milky-Way-like potential may vary between 
approximately 0.5ft and 2ft. Fortunately, the flat rotation 
curve of the Milky Way implies that there is a wide range 
of galactocentric distances within the disk where k is about 

k~ 1.4ft, (11) 

(cf. Fig. 3 in lJust et al1l2009h . 

As mentioned above, this theory only holds for nearly 
circular orbits. In the case of arbitrary orbits in complicated 
potentials these approximations break down and no simple 
solutions to the equations of motion can be found. 

In the following analysis the tidal radius is always de- 
termined by combining equation and [TU1 i.e. 

GM 



3 _ 

L ~ ft 2 - d 2 $>/dR 2 ' 
The angular velocity is furthermore calculated using 



ft = 



\R x V\ 
R 2 ' 



(12) 



(13) 



where R is the position of the cluster with respect to the 
galactic centre and V its velocity. 

The basic quantities which we investigate (bound mass, 
mean stellar mass, etc.) will be calculated for all stars within 
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this theoretical radius, and the distance of the detected over- 
densities relative to the density centre compared to it. This 
distance, yc, should only be related to the tidal radius if the 
overdensities are of epicyclic origin. 

2.5 yc in the case of time-dependent tidal fields 

There is no theory for the distance of the epicyclic overden- 
sities, yc, in the case of time-dependent tidal fields yet. But 
we can estimate the distance between the cluster and the 
epicyclic overdensities, even when the system gets period- 
ically accelerated and decelerated on an eccentric orbit or 
through a galactic disk, if we make further assumptions: 

(i) Firstly, we have to consider the fact that the acceler- 
ation of the cluster and of the stars in the tails will change 
with time. Consequently, the distance between any two 
points moving along the orbit is not constant but changes 
according to the orbital phase. What stays constant due to 
angular momentum conservation though, is the difference 
between the times at which the two points pass a certain 
point of the orbit. We can therefore convert the distance yc 
into a conserved time difference At. For constant tidal fields 
yc can be expressed as 



yc = VAt, 



(14) 



where V is the magnitude of velocity of the cluster on its 
circular orbit about the galaxy and At is the time difference 
at which the cluster and the overdensity pass a certain point 
of the orbit. 

(ii) For eccentric orbits or orbits including disk shocks 
the cluster and its tails get periodically stretched and com- 
pressed as the velocity of the cluster and the tails changes 
during a period. The distance between two points which are 
separated by a given amount of time, At, can be written in 
a Taylor expansion as 



yc(t) = V(t)At + -A(t)At 2 + 0(At 3 ), 



(15) 



where A is the acceleration at the given point. Since V and A 
are well known for our models, the separation of two points 
which would be at a constant distance in a constant tidal 
field, e.g. the cluster and the first-order overdensity, can be 
calculated for a time-variable tidal field. 

(iii) Furthermore, we shall assume that the time varia- 
tions of the tidal field are sufficiently fast that the cluster 
cannot adapt to the changing environment, but rather be- 
haves as if it experiences a single mean tidal field along its 
orbit. For this purpose we shall average the galactocentric 
distance, R, and the orbital velocity, V, and use these quan- 
tities to calculate an average yc using eq.[8]and assuming a 
circular orbit at the given mean galactocentric distance. 

The results of these assumptions will be checked in 
Sec. 13.41 but the following argument explains why assump- 
tion (iii) is reasonable. The length of the epicycle, which 
escaping stars follow after they leave the cluster vicinity, 
depends on the offset of the star from the cluster centr e per- 
pendicular to the orbit (see KMH or I Just et aL I l2009l for a 
more detailed discussion). In a constant tidal field this off- 
set is the tidal radius. But for time- variable tidal fields, the 
stars do not necessarily have to escape from the tidal ra- 
dius, as this radius might change faster than the cluster can 
adapt to the changing tidal conditions. In fact, stars escape 



from what we might refer to as the "edge" of the cluster, 
whose radius changes little during one orbit. Indeed, recent 
studies show that stars escape preferentially from a radius 
approximately equal to the apogalactic tidal radius of the 
cluster; furthermore they have low velocities, making the for- 
mation of epicyclic overdensities likely, even in the case of 
non-circular orbits (Kiipper, Kroupa, Baumgardt & Heggie, 
in preparation). 

(iv) Finally, by applying eq.Q3]we can estimate a time by 
which the cluster and the overdensity should be separated 
on their orbit. With this approximate value of At we have 
all the ingredients to predict the behaviour of yc with time 
during one orbital period using eq. 1151 



2.6 Mean drift velocity 

The time escaping stars need to complete one epicycle, e.g. 
the time that passes from the time a star leaves the cluster 
until it reaches the first overdensity, is given by 



tc = 2it/k, 



(16) 



and in KMH we also showed that the distance between the 
overdensities is given by equation[8] Thus, we can determine 
a mean drift velocity of the stars within the tidal tails with 
respect to their guiding centre by using 



vc 



yc_ 

tc 



±2n l 



4» 2 

*2 



XL, 



(17) 



where the minus (plus) sign holds for the leading (trailing) 
tail. Hence, the mean drift velocity depends on the tidal-field 
properties through k and xl- This equation is only valid for 
circular orbits, though. For circular orbits in the galactic 
disk of the Milky Way the mean drift rate is about 

v c ~±2flx L = ±{AGMQ.) 1/3 , (18) 

where we used equations QIJ] and QT] For a constant mass- 
loss rate the density within the tails is therefore increasing 
with decreasing mass as the angular velocity of the cluster is 
constant whereas the tidal radius decreases monotonically. 
Hence, the stars propagate more slowly along the tails the 
less massive the cluster is. 

In time-dependent tidal fields the mean tidal conditions 
used in the formalism of Sec. 12.51 may be applied to get a 
rough estimate of the mean drift velocity. But due to the 
compression and stretching of the tails throughout the orbit 
a drift velocity in a literal sense is not meaningful. 

The mean drift velocity, vc, which is calculated in the 
rotating reference frame of the leading or trailing tail, can 
be converted into the average orbital velocity of Sec. 12.31 
which is measured in_the reference frame of the cluster, by 
making a correction, V , for the difference in orbital velocity 
of the guiding centre of the corresponding tail to the orbital 
velocity of the cluster. This correction emanates from the 
left part of the right-hand side of eq. 4 in KMH, which gives 
the offset of the guiding centre of the tails from the cluster 
centre, multiplied by the angular velocity fl, and reads 



V 



+o 4n2 

±S2 7T-XL, 

K 



(19) 



where the plus (minus) sign holds for the leading (trailing) 
tail. The average orbital velocity in the tails with respect to 
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Table 1. Overview of all computed models. In the first column 
the section is given where the particular model is discussed. R apo 
gives the apocentre distance from the galactic centre and R per i 
the pericentre distance, incl the inclination of the orbit with 
respect to the galactic disk and e the eccentricity of the orbit 
(eq. 126 I t - All clusters have initially 65536 stars and a mass of 
about 20000 Mq. The models with inclination of 90 deg to the 
disk always cross the disk at R '■ the eccentricity of these orbits 
is simply due to the asymmetric galactic potential. 



Sec. 


Rapo [kpc] 


Rperi \kpc\ 


incl [deg] 


e 


EI 


8.50 


8.50 


_a 


0.00 


|q q| 
lo.ol 


O. ou 


O. ou 


u 


u.uu 


13.41 




5.10 





0.25 






2.83 





0.50 






1.21 





0.75 


l3~5l 


9.27 


8.50 


90 


0.04 




4.25 


4.25 





0.00 






2.55 





0.25 






1.42 





0.50 






0.61 





0.75 




5.24 


4.25 


90 


0.10 




12.75 


12.75 





0.00 






7.65 





0.25 






4.25 





0.50 






1.82 





0.75 




13.27 


12.75 


90 


0.00 




17.00 


17.00 





0.00 






10.20 





0.25 






5.67 





0.50 






2.43 





0.75 




17.44 


17.00 


90 


0.01 



a Cluster is moving in a point-mass galactic potential. 

the cluster, tJ, is then given by 

~ / 4fi 2 \ 4fi 2 
v = vc + V = ±2tl[l T )x L ±n—XL (20) 

= ±n(^--2^x L . (21) 

This equation gives the mean value of the velocity profile 
along the tidal tails which will be measured in the following 
investigation using eq. [5] 



3 NUMERICAL INVESTIGATION 

The aim of this section is to trace the behaviour of tail for- 
mation from the ideal case of KMH to the complex case of a 
realistic star cluster. Therefore, we start with a star cluster 
like the one of KMH and, as a first step, add a mass spec- 
trum to the cluster stars (Sec. I3.2|l to see if it affects the 
formation of epicyclic overdensities. Thereafter, we change 
the point-mass galactic potential to a Milky- Way potential 
with bulge, disk and halo (Sec. I3.3|l . This will then be our 
reference model for the rest of this paper, i.e. we shall com- 
pare each model with this one. Afterwards, we study the 
influence of different orbital types on the formation of sub- 
structure. First, the orbit is changed from circular to eccen- 



tric (Sec. 13. 4p and then a set of computations is made with 
circular orbits but perpendicular to the disk (Sec. I3.5|l to 
get an idea of the importance of disk shocks on cluster evo- 
lution. Additionally, all these kinds of orbits are computed 
at galactocentric distances of 4250, 8500, 12750 and 17000 
pc. First the clusters at an initial galactocentric radius of 
8.5 kpc are discussed, but the effect of changing this radius 
is studied separately in Sec. 13.61 An overview of the models 
can be found in Table [T] 

3.1 The setup 

The initial clusters are chosen such that they can represent 
the Milky- Way globular cluster Palomar 5 (which is itself 
about 10 Gyr old) during the last 3 Gyr of evolution, i.e. 
rather 'fluffy' wit h a ratio of half-mass radius to tidal ra- 
dius of about 0.2 (lHarrisI 1996). We do this because Pal 5 is 
the only globular cluster with prominent tidal tails in which 
significant s ubstr ucture has been found (|Odenkirchen et al.l 
l200ll |2002| . l2003h . Furthermore, there are hardly any clus- 
ters present in the Milky Way which are more vulnerable 
to tidal influences than Pal 5. Most globular clusters in the 
Milky Way are much more compact and therefore can resist 
pericentre passages or disk shocks much better. Hence, the 
complete set gives us an idea of the maximum importance 
of a wide range of realistic tidal influences on star clusters 
and their tidal tails. 

All computations were perfo rmed wi t h th e colli - 
sional JV-body code NBODY4 (|Aarsethl 1 19991 . I2003T I 
on the GRAPE-6A supe r comp uters at the AlfA 
jFukushige. Makino fc Kawail l2005h . The internal ini- 
tial conditions of the chosen clusters are all similar: the 
clusters contain 65536 stars drawn fr om a canonica l , two - 
part power law initial mass function l|Kroupa|[200ll . 120081 ) 
ranging from 0.1 to 1.2 Mq; this results in a total cluster 
mass of about 20000 Mq, which is c omparable to t h e mas s 
of Pal 5 3 Gyr ago as estimated by iDehnen et al.l l|2004h . 
The stars follow a Plummer density profile with a half-mass 
radius chosen such that Rh/Rude = 0.2 for a circular orbit, 
i.e. all clusters at a given apogalactic distance use the same 
cluster setup to ensure comparability among each other. 

The clusters are set up as isolated clusters and then 
put into the tidal field, such that there is a tide-dependent 
unbound fraction of stars at the very beginning, which is in 
most cases negligible and mostly results in a slightly higher 
mass-loss rate at the beginning of the computations. 

The clusters do not contain any primordial binaries be- 
cause binaries would significantly enhance the CPU demand 
without adding to the problem at hand. Furthermore, stellar 
evolution is neglected to focus on dynamical effects. After 
setup the calculations were carried out for 4 Gyr if possible. 

Only one cluster was computed for a point-mass galac- 
tic potential; for all other clusters a Milky- Way potential 
consisting of bulg e , disk and spherical halo as described in 
lAllen fc Santillanl l|l99ll ) was used. For the latter potential, 
the equation for x l (eq. 1 10 p cannot be used unrestrictedly, 
because it only holds for spherical symmetric potentials in 
which the total angular momentum is naturally conserved 
or for orbits lying in the plane of symmetry, i.e. for orbits in 
the disk, for which the component of the angular momen- 
tum perpendicular to the disk is conserved. Nevertheless, 
even for the orbits which do not lie in the plane of symme- 
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Figure 2. Number distribution of stars along the tidal tails in 
bins of 25 pc (upper panel) and corresponding mean velocities 
of the stars within the bins (with respect to the orbital velocity 
of the cluster, eq. [5} for the cluster in the point-mass galactic 
potential (discussed in Sec. I3.2C The snapshot was taken when 
the cluster had a mass of about 16300 Mq at t = 2.0 Gyr, which 
means a theoretical tidal radius of roughly = 33 pc and a 
predicted distance of the first-order overdensity of 1250 pc. Error 
bars in the upper panel are Poisson errors and in the lower panel 
they show the velocity dispersion within the corresponding bin 
(eq.[6j. 



try the equation is valid to a sufficient degree of accuracy 
because the variation of the total angular momentum is only 
of order 10%. 



3.2 Point-mass galaxy, circular orbit 

The step from a cluster consisting of single-mass stars to a 
multi-mass cluster is necessary, although no deviation from 
the theoretical predictions of KMH is expected, since the so- 
lutions to the equations of motion of stars leaving a cluster 
through its Lagrange points (equations 4-6 therein) do not 
depend on the masses of the stars. The only thing to check is 
whether the evaporative escape conditions are still fulfilled. 
But, as discussed in KMH, in constant tidal fields this con- 
dition is better fulfilled the more stars there are within the 
cluster, and here we are dealing w ith 64k sta r s com pared to 
N = 1000 in KMH. Furthermore, lJust et all (|2009h already 
proved the formation of epicyclic over- and underdensities 
for this kind of cluster. Nevertheless, we perform this exper- 
iment to confirm our methodology. 

Fig. [2] shows a snapshot of the system at t = 2.0 Gyr 
when the cluster has a bound mass of about 16300 Mq , i.e. a 
tidal radius of 33 pc. The cluster clearly shows overdensities 
at the predicted positions of about yc = ±12ttxl m 1250 pc 
and even shows signs of the second- and third-order overden- 
sities at ±2t/c and ±3t/c respectively (which are not visible 
in the plot because the plotted range of y is restricted to 
facilitate comparison with similar plots for other clusters). 
Also the orbital velocity with respect to the cluster in the 
lower panel of Fig. [2] clearly shows the predicted signal of 
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Figure 3. Mass evolution of the clusters at an initial galacto- 
centric distance of 8.5 kpc. Here M is the mass inside radius 
xl , determined using equation HOI The reference cluster is the 
most unperturbed and therefore survives for the longest time. 
The point-mass model compares well with the reference model 
since they have the same concentration and are both not subject 
to tidal variations. The cluster with the inclined orbit (Sec. I3.5H 
shows the slight influence of disk shocks, which are rather neg- 
ligible at a galactocentric radius of 8.5 kpc for a cluster of the 
given concentration. On the contrary, eccentric orbits (Sec. I3.4H 
significantly decrease a cluster's lifetime. Also visible from the 
figure is the fact that for all models mass-loss takes place rather 
smoothly and does not happen in steps. Moreover, the effect of 
the pericentre dips of the tidal radius on the bound mass can be 
seen, but it shows that most stars get recaptured by apocentre. 



periodic acceleration and deceleration. Moreover, the aver- 
age orbital velocity in this figure is in good agreement with 
the predictions of ea. l20l yielding an average orbital velocity 
of about v = ±1.7kms -1 with respect to the cluster for the 
values of Q and xl given in the caption of Fig. [2] 

From this figure we can also see that most escaping stars 
have a velocity which is larger than the escape velocity, by 
the following argument. Eq. [8] gives the distance of the first 
epicyclic loop for a star which passes the Lagrange point 
with vanishing velocity; this distance is about 1250 pc in 
Fig. [2] If a star has an excess velocity, v, at the moment of 
es cape then the le ngth of the epicycle increases (see eq. 22 
in lJust et ai]l2009T l as 



yc = ±- 



4ttD. 



4tf 

*2 



XL + 



V 

2Q 



(22) 



In Fig. [2] the first epicyclic maximum is at about 1500 pc. 
This is consistent with eq.[22]if we assume that the average 
excess velocity is about 0.34 km/s. 

3.3 Milky- Way potential, circular orbit 

The cluster described in this section will be the reference 
for all following clusters, since all of them are moving in the 
same galactic potential but on more complex orbits. Here 
we discuss some basic evolutionary parameters of this clus- 
ter (i.e. bound mass, Lagrange radii and mean stellar mass) 
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Figure 7. Time evolution of the number distribution of stars along the tidal tails in bins of 25 pc for the reference cluster (discussed 
in Sec. I3,3j l, The first- and second-order overdensities can be seen at about ±400 pc and ±800 pc, respectively, decreasing slightly in 
distance as the cluster loses mass and the tidal radius decreases. The color bar on the right shows the color coding of the number of stars 
per bin. 
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Figure 4. Evolution of the 2%, 10%, 50% (R h ) and 90% Lagrange 
radii of the reference cluster. Additionally, the tidal radius (ij) 
and the core radius (Rc) are shown. The former shows the ex- 
ternal/tidal evolution of the cluster which, in this case, is just 
affected by the mass-loss of the cluster itself, while the latter re- 
flects the internal dynamical evolution of the cluster. By 4 Gyr 
the cluster is beginning to go into core collapse. 



and only point out significant differences for the other clus- 
ters in the following sections. The aim is to understand the 
evolution of the reference cluster in detail before we look at 
the tidal tails and the overdensities therein. 
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Figure 5. The average mass of stars, M/N, versus bound mass, 
M, for the clusters at 8.5 kpc. The preferential loss of low-mass 
stars is very pronounced in all clusters. Since preferential loss is a 
phenomenon of two-body relaxation, the most unperturbed clus- 
ters, i.e. the reference cluster and the cluster in the point-mass 
potential, show the earliest increase in M/N with decreasing M. 
As all clusters, except the point-mass case, have the same initial 
relaxation time, the start of mass segregation gives a quantifica- 
tion of the importance of two-body relaxation compared to the 
tidally induced mass loss. Computations were halted after 4 Gyr 
so not all clusters dissolve completely. 
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Figure 6. Number distribution of stars along the tidal tails in 
bins of 25 pc (upper panel) and corresponding mean velocities of 
the stars within the bins (with respect to the orbital velocity of 
the cluster, eq. O for the reference cluster (discussed in Sec. 13 - 3 1 ) - 
The snapshot was taken when the cluster had a mass of about 
16400 Mq at t = 2.0 Gyr, which means a theoretical tidal radius 
of roughly 37 pc and a distance of the first-order overdensity 
of about 350 pc (equation [8} . Error bars in the upper panel are 
Poisson errors and in the lower panel show the velocity dispersion 
within the corresponding bin (eq. [6}. 



3.3.1 Internal evolution 

Putting the cluster of Section r3.2l into a Milky-Way potential 
will not significantly change its evolution as long as it still 
moves on a circular orbit in the disk and with the same 
angular velocity. The only change is the steepness of the 
potential which, in the formalism described above, alters the 
epicyclic frequency, k 2 , and therefore changes the value of 
the tidal radius. Since the Milky- Way potential in the disk 
falls off less steeply, d 2 $/dR 2 is smaller and the tidal radius 
therefore larger. 

This cluster is not perturbed by a time-dependent tidal 
field. Hence, as no external energy comes in, its lifetime is 
longest (Fig [3} . During the first 500 Myr there is a small 
burst of primordial escapers. Thereafter, the mass- loss rate 
is about constant. After four Gyr the cluster still has about 
75% of its initial mass and is only just beginning to go into 
core collapse (Fig. Q. The core radius, Rc, in this figure 
shows the internal dynamical evolution of the reference clus- 
ter. It is calculated using the algorithm (cf. lCasertano fc Hutl 

GUI)) 



Rc = 



N 2 2 



Ere 2 
i Pi 



(23) 



where ri is the distance of the ith star from the density 
centre and pi is the mass density around this star, computed 
using the distance to its fifth closest neighbour and the mass 
within this radius. 

A good indicator for the formation of epicyclic overden- 
sities is ongoing two-body relaxation in the cluster, as escape 
due to two-body relaxation ensures evaporative escape con- 



ditions (Sec. 2. 2). In Fig. [5] the mean mass of bound stars 
versus bound mass is shown. Since a cluster is expected to 
lose preferentially low-mass stars, as those are most proba- 
bly accelerated to velocities above the escape velocity, and 
this effect is due to energy equipartition (via two-body re- 
laxation), Fig. Ogives a measure of the effectiveness of two- 
body relaxation for the escape of stars. The preferential loss 
is clearly visible in all clusters. We can see that the reference 
cluster, as it is the most unperturbed one, shows the earliest 
increase in mean mass with decreasing bound mass. 



3.3.2 Evolution of the tidal tails 

As we have seen that two-body relaxation plays a major role 
in the dissolution process, and escape due to two-body re- 
laxation leads to epicyclic overdensities, the reference cluster 
also shows the expected density enhancements at the pre- 
dicted positions (Fig. [6}. As mentioned in Sec. 12.41 we can 
use k, = 1.4fi for this orbit in the Milky- Way potential, 
which gives 



yc 



GM 
2?F' 



(24) 
(25) 



Hence, the epicycles of the escaping stars are about four 
times shorter than in the case of a point-mass galaxy while 
the tidal radius is almost 10% larger. Still, the epicyclic ap- 
proximation is correct, as can be seen in Fig. |S] where the 
binned stars along the tidal tails are shown at 2.0 Gyr when 
the cluster has a bound mass of about 16400 Mq. The corre- 
sponding tidal radius and the predicted position of the first 
overdensity, yc, are given in the figure caption. As in Fig. [2] 
the first-order epicyclic maximum is at a somewhat larger 
distance than the predicted yc- The difference is about 50- 
100 pc which, using eg. 1221 corresponds to an average excess 
velocity of about 0.27-0.55 km/s. 

The lower panel of Fig. HJ] also shows the predicted ve- 
locity variations. Moreover, the average orbital velocity with 
respect to the cluster is again in good agreement with the 
predictions of eq. 1201 which yields an average orbital velocity 
of about v — kms" 1 for the given parameters. This average 
value is intuitively expected, since the rotation curve of the 
assumed galactic potential is almost flat at this galactocen- 
tric radius. The role of the fiat rotation curve can be nicely 
seen by comparing Fig. [5] with Fig. [5] where larger velocity 
differences between the leading and trailing tidal tails are 
seen. 

When comparing Fig.|S]with Fig. [5] which show similar 
clusters at the same dynamical state, i.e. at a bound mass 
of about 16000 Mq, but in a different tidal field, it should 
be noticed that, besides the fact that yc is approximately 
four times larger in Fig. [2] the overdensities are also far 
more extended in Fig. [2] also by a factor of about four. 
Furthermore, the peak density outside the cluster is larger 
by a factor of two in the case of the more realistic reference 
cluster. This is due to the mean drift velocity of the stars 
along the tidal tails (eq. I17p . Comparing vc for the two 
clusters we see that, while the escaping stars of the cluster 
in the point-mass potential have a mean drift velocity of 
about 5 kms -1 , the tail stars of the reference cluster move 
with only about 2 kms -1 along the tail, i.e. the tail is denser 
on average because the mass-loss rate of both clusters is 
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Table 2. Parameters for the formalism described in Sec. 12.51 for 
the models at 8.5 kpc distance. The first two columns are taken 
from Table [T1 whereas R gives the mean galactocentric radius, V 
is the mean orbital velocity of the cluster and At is the estimate 
of the time difference between the cluster and the first-order over- 
density reaching the same phase of the orbit on their orbit about 
the galaxy as derived in Section [23] 
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approximately the same (see Fig. [3J). This implies that the 
tails of the reference cluster, and especially its overdensities, 
would be much easier to observe. 

Fig. [7] shows a time series of the number distribution 
along the tidal tails. In this figure the cluster is omitted to 
allow a larger dynamical range in the representation of the 
number density of the tails. The first- and second-order over- 
densities can be seen with distances decreasing slightly with 
time as the mass, and hence the tidal radius, of the system 
monotonically decreases. There is also an initial transient 
overdensity, which propagates rapidly away from the cluster 
at small times. This is due the initial presence of "primordial 
escapers" (Sec. 3.1), and is an artefact of the initial condi- 
tions. But it is an interesting feature, as it illustrates how 
very differently the overdensities evolve if they are due to a 
single pulse of mass loss. 

3.4 Milky- Way potential, eccentric orbits 

We study the influence of pericentre passages on the forma- 
tion of tidal tails by using eccentric orbits, but remain within 
the disk in order to avoid disk shocks. The eccentricities, e, 
of the chosen orbits are 0.25, 0.5 and 0.75, where 

Rapo Rperi 



Rapo ~t~ Rpe 



(26) 



The clusters start at a radius of 8.5 kpc and have an ini- 
tial velocity which is reduced by a corresponding amount 
compared to the circular velocity, and hence 8.5 kpc is their 
apocentre distance. The corresponding values of R per i can 
be found in Table [1] 

3.4-1 Internal evolution 

The clusters computed in this way dissolve faster than the 
reference cluster (Fig. [3}, which is expected, because the 
more eccentric the orbit the stro nger is the average tidal 
field (Bau mgardt fc Makinol 12003), Furthermore, for an ec- 
centricity of 0.25 there are already small periodic wobbles, 
which are visible in Fig.|3]and which grow with increasing e. 
These are due to the periodically changing value of the tidal 
radius xl, which is smallest at pericentre, and this sets loose 
a fraction of stars, of which most are recaptured by apoc- 
entre. As discussed in Section 12.1.31 due to this recapture, 
stars are not lost in a pulse during a pericentre passage: the 
mass curve rather follows a smooth line (if the pericentre 
dips are neglected), and not a series of steps. 
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Figure 9. Distance of the first-order epicyclic overdensity, yc, 
estimated with the formalism described in Sec. I2.5l for the eccen- 
tric orbits at 8.5 kpc. The corresponding values of R, V and At 
can be found in Table [2] The values for the velocity V and the 
acceleration A on the cluster were extracted from our computa- 
tions in time steps of 10 Myr. For comparison the value of the 
reference cluster is also shown. On eccentric orbits the distance 
yc starts oscillating about a mean value, where the amplitude of 
the oscillations increases with increasing eccentricity. The mass 
of the cluster in this plot is fixed to the initial mass of 20000 
Mq. If mass loss was taken into account the mean value as well 
as the amplitude of the oscillations would decrease with time. 
The predictions made within this figure can be compared with 
the Af-body results in Fig. Upl and Fig. 1111 (We can also see from 
this figure that, unlike in a Keplerian potential, the period of 
the orbits decreases with increasing eccentricity. This further in- 
creases the effect of higher eccentricities on the internal cluster 
evolution.) 



This behaviour, in which we follow the mass within the 
tidal radius calculated using eg. 1101 i s very different from 
the behaviour found by iDehnen et al.l (120041 ) (their fig. 4), 
who included only stars lying within the initial tidal ra- 
dius for the computation of internal quantities of the clus- 
ter, such as the total mass. Their approach suggests that 
stars are lost in pulses during pericentre passages or disk 
shocks, which is in strong contrast with our finding. Differ- 
ences in the initial conditions, orbital parameters, galactic 
model and even computational method all have a role in this 
comparison. We note, however, that recent results (Kiipper, 
Kroupa, Baumgardt & Heggie, in preparation) show that 
for computations of internal quantities it is more appropri- 
ate to use the apogalactic tidal radius or the mean tidal 
radius mentioned in Sec. 12.51 

While the strength of these dips in the mass curve do 
not reflect permanent escape, they do indicate what frac- 
tion of stars is seriously affected by a pericentre passage or 
a disk shock, since it shows how many stars are at large 
radii and get temporarily unbound. (Those stars which are 
most affected by a pericentre passage or a disk shock are 
those at large radii (cf. equations [3] and [4| ) . But, as we can 
see in these dips, not all stars that have been once outside 
the tidal radius are lost from the cluster. As mentioned in 
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Figure 8. Time evolution of the mass shells (Lagrange radii) of the cluster with an eccentricity of 0.5 and an apocentre distance of 8.5 
kpc. The lines show the radii containing 2, 5, 10, 20, 30, 40, 50, 60, 70, 80 and 90 percent of the initial cluster mass (about 20000 M©). 
The periodically changing line, which is uppermost at t = 0, is the tidal radius. The periodic expansion and compression of the whole 
cluster through the pericentre passages can be observed in all shells. Only the innermost shells show a progressive contraction while the 
other shells show an accelerating expansion. If a shell gets unbound during a pericentre passage, the shell will quickly expand and move 
away from the cluster. If the shell can get recaptured by the tidal radius, the expansion will be slowed down or even reversed. 



Sec. 12.1.31 the loss of these outlying stars cannot happen 
instantaneously, essentially because their orbital time about 
the cluster is on the order of the orbital time of the cluster 
about the galactic centre. 

The finite breadth of the dips (e.g. for e = 0.25) also 
tells us that not all eccentricities lead to tidal shocks in a 
classical sense. Therefore, we have tried to avoid using the 
term shock in connection with eccentric orbits, and speak 
of the more general pericentre passages, which implies a 
broader range of time scales. 

The behaviour we are describing can be observed in 
Fig. [S] where the time evolution of mass shells of the clus- 
ter with an eccentricity of 0.5 and an apocentre distance 
of 8.5 kpc is shown. The figure demonstrates the periodic 
expansion and compression of the whole cluster, but the 
amplitude is much smaller than the change in tidal radius. 
Superimposed on this externally induced oscillation is the 
internal dynamical evolution of the cluster. While the in- 
nermost shells are contracting due to mass segregation, the 
remaining shells are expanding. 

If a shell gets unbound during a pericentre passage it 
will expand thereafter, where the amount of expansion de- 
pends on the time the shell spends outside the tidal radius. 



If this time is short (e.g. the curve for 80% in Fig. [H] at 
about 750 Myr, the expansion will be slow enough that the 
shell can be recaptured by the growing tidal radius when the 
cluster moves from peri- to apocentre. A recapture will de- 
celerate the expansion or even reverse it. When no recapture 
is possible, the shell will quickly move away from the cluster 
vicinity within a fraction of the cluster orbital period. Here 
again it can be seen that escape happens throughout the 
whole orbit and is not heavily concentrated at perigalacti- 
con. 

But how does the eccentricity of the orbit affect two- 
body relaxation in the cluster and what is the effect on the 
escape conditions? Fig. [5] shows that preferential loss of low- 
mass stars is increasingly suppressed compared to the refer- 
ence cluster with increasing eccentricity, since the fraction 
of lost stars with arbitrary masses increases. Nevertheless, 
all clusters show an increasing mean stellar mass, such that 
two-body relaxation is evidently at work. Provided that the 
escape conditions for a majority of stars are still evapora- 
tive, and the eccentricity of the orbits does not induce too 
large a large scatter in escape conditions, we may still ex- 
pect that epicyclic overdensities will arise. We now examine 
the tidal tails to check this expectation. 
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Figure 10. As Fig. [7] time evolution of the number distribution of stars along the tidal tails in bins of 25 pc for the cluster with 
apogalactic distance of 8.5 kpc and orbital eccentricity of 0.25 (discussed in Sec. 13.3b . The first- and second-order overdensities can be 
seen, periodically changing in distance but in mean following the trend of the reference cluster (cf. Fig. 0. The computations agree with 
our predictions made in Sec. [23] fairly well, as can be seen when compared with Fig. [9] 
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Figure 11. As Fig. [7] time evolution of the number distribution of stars along the tidal tails in bins of 25 pc for the cluster with 
apogalactic distance of 8.5 kpc and orbital eccentricity of 0.5 (discussed in Sec. |3.3t . The first-order overdensities can be seen especially 
in perigalacticon as a faint blue structure, following our predictions made in Sec. 12.51 fairly well (cf. Fig. [5}. The most prominent 
overdensity at early times, however, is due to primordial escapers. It can be followed as a red structure moving slowly away from the 
cluster. 
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Figure 12. Number distribution of stars along the tidal tails in bins of 25 pc (upper panels) and corresponding mean velocities of 
the stars within the bins (with respect to the orbital velocity of the cluster, eq. [5j for the cluster with an eccentricity of 0.5 and an 
apogalactic distance of 12.75 kpc (discussed in Sec. 13.61 1. The twelve consecutive snapshots were taken when the cluster had a mass of 
about 14000 M starting at t = 970 Myr in time steps of 10 Myr. A full revolution takes about 170 Myr. So, the first panel is 20 Myr 
after the cluster has been at perigalacticon where the tails where maximally stretched, in the middle panel of the third row (t = 1040 
Myr) the cluster is in apogalacticon and maximally compressed. In the right panel of the forth row the cluster is 40 Myr away from 
perigalacticon. In this time series the difference in average number of stars per bin between peri- and apogalacticon is observable which 
exceeds a factor of three (compare, for example, the left panels in the first and third row in the vicinity of the cluster). Furthermore, the 
asymmetric compression can be seen, wandering from positive to negative y- values (indicated by the arrow), as an asymmetric density 
distribution when the leading and the trailing tail are carefully compared. The part of the tails which is compressed most shows a number 
density which is about twice as high as the same part in the opposite tail. The asymmetric acceleration and deceleration of the tail stars 
can also be seen in the velocity distributions, where the velocity dispersion is large in parts of the tails which are broad and small in 
parts which are thin. 
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3.4-2 Evolution of the tidal tails 

The tidal tails of clusters on eccentric orbits behave dif- 
ferently from those of clusters on circular orbits in several 
ways. First of all, the coordinate system in which the cluster 
is at rest is no longer uniformly accelerated, and therefore 
the epicycle approximation cannot be applied without mod- 
ification. For this reason, the solutions of the equations of 
motion of KMH (therein eq. 4-6) do not hold any more. 
Hence, there is no simple relation between the tidal radius 
and the distance of the density maxima unless we make fur- 
ther assumptions (Sec. [275)) . Nevertheless, the equation for 
the tidal radius (eq.llO[) is still valid since the orbits lie in the 
plane of symmetry of the potential and therefore the angu- 
lar momentum is conserved. But the value of d 2 &/dR 2 now 
changes with time as the galactocentric distance oscillates, 
and the angular velocity, of course, also varies with time. 
This gives a variation of the tidal radius (between peri- and 
apocentre) of order 90% of the apogalactic value for e = 0.75 
for the clusters with an apogalactic distance of 8.5 kpc. 

Epicyclic overdensities 

As we expected from theoretical considerations above, we 
also find epicyclic overdensities in the tidal tails of clusters 
on eccentric orbits. Indeed, all clusters in the range of eccen- 
tricities we studied show these overdensities. The epicyclic 
origin of the overdensities we have found can be seen by 
looking at time series of the number distribution along the 
tidal tails, as the distance yc should be related to the tidal 
radius. 

Since the tidal radius is large at apocentre and small 
at pericentre, the density maxima are expected at a larger 
distance at apogalacticon and vice versa. But the opposite 
is the case, which is due to the non-linearly accelerated rest 
frame: from apo- to perigalacticon the cluster and its tails 
get accelerated and stretched, while from peri- to apocentre 
the tails g et compressed as the the whole system is decel- 
erated (cf. iPiatek fc Prvorlll995l ). Hence, at pericentre the 
cluster is fastest and the distance between the maxima is 
largest, whereas the maxima and the whole tails get most 
compressed at apogalacticon - which is also the explana- 
tion why the tidal tails of Palomar 5 are so clearly visible 
since it is near the apoce ntre of its orbit about the Galaxy 
jOdenkirchen et al.ll^Ool ). 

If we apply the assumptions of Sec. 12.51 we can predict 
the positions of the overdensities even for the eccentric or- 
bits. Therefore we average the galactocentric radius and the 
orbital velocity of the cluster for each model in Table [2] and 
estimate a time difference, At, by which the cluster and the 
first-order overdensity should be separated in reaching the 
same phase of the orbit. The variations of yc for the differ- 
ent orbital types at a fixed mass of 20000 Mq are shown in 
Fig- HI for several revolutions. 

Obviously, the distance yc now strongly depends on 
the position of the cluster on its orbit; from apo- to pericen- 
tre it varies increasingly with increasing eccentricity. This 
periodic change of yc is also observable in Fig. [10] where 
the change of the number distribution along the tidal tails 
with time is shown for the cluster at an apogalactic distance 
of 8.5 kpc and an eccentricity of 0.25. While the reference 
cluster (Fig. [7} shows the predicted monotonic decrease of 



yc, the clusters on eccentric orbits have an orbital variation 
imprinted on this monotonic decrease. Our estimate from 
Sec. 12.51 shows good agreement with the simulations; both 
vary by a factor of 1.7 about a mean value which is lower 
than the corresponding value of the reference cluster. The 
same holds for the cluster with an orbital eccentricity of 0.5, 
shown in Fig. 1111 The epicyclic overdensities in this figure 
are quite faint compared to the initial mass loss, which is 
due to primordial escapers. Nevertheless, the location of the 
first-order epicyclic overdensity is in good agreement with 
the predictions made using the formalism of Sec. I2,5F1 . sup- 
porting the assumptions made for this purpose. 

While yc in Fig. varies for an eccentricity of 0.25 by a 
factor of about 1.5 between maximal and minimal distance 
from the cluster centre, and for an eccentricity of 0.5 by a 
factor of 2.7, the amplitude of variation for e = 0.75 is even 
a factor of 7.5. Moreover, for the latter eccentricity the com- 
pression at apogalacticon is so strong that the overdensities 
are literally pushed back into the cluster such that they over- 
lap with each other and the cluster and nearly disappear in 
our number distribution plots. 

Asymmetry overdensities 

The non-uniform acceleration along the orbit also leads to 
an asymmetry in extended tidal tails, since in the case when 
the cluster is exactly at apocentre, for example, the leading 
tail is already accelerating towards the galactic centre while 
the trailing tail is still slowing down. The difference in ac- 
celeration results in an asymmetry in compression between 
the leading and the trailing tail, and thus in an asymmetric 
stellar density within the two tails. The asymmetry increases 
with a larger eccentricity and is also visible in the velocity 
distribution along the tidal tails. 

Nevertheless, the epicyclic overdensities stay visible 
throughout the whole orbital period. We demonstrate this 
for the cluster with an orbital eccentricity of 0.5 but at an 
apogalactic distance of 12.75 kpc. We do this because we 
generated outputs of NBODY4 every 10 Myr throughout 
our computations, and for clusters at smaller apogalactic 
distances the time resolution is too poor for a convincing 
demonstration. For this cluster with e = 0.5 the orbital time 
is about 170 Myr, we therefore show 12 consecutive snap- 
shots starting at 970 Myr in steps of 10 Myr in Fig. [12] At 
that time it has a mass of about 14000 Mq and is first close 
to pericentre, then at 1040 Myr in apocentre and in the last 
panel, at 1080 Myr, on its way to pericentre again. 

Not only is the compression significant, increasing the 
stellar density within the tails by, at least, a factor of three 
at apogalacticon compared to perigalacticon, but also the 
asymmetry is clearly visible and can be followed throughout 
the snapshots wandering from right to left (indicated by a 
black arrow) and enhancing the number density by a factor 
of about two compared to the same part in the opposite tail. 

2 For the cluster with an orbital eccentricity of 0.75 the time 
resolution of our output data is too poor compared with the or- 
bital time. Moreover, the average density within the tails differs 
so much between apo- and perigalacticon such that it is hard 
to convincingly show the epicyclic overdensities in this kind of 
representation. But this is just a problem of representation. 



Tidal Tails of Star Clusters 15 



o 

Q. 



O 



1000 



800 



600 



400 



200 




400 600 
t [Myr] 



1000 



Figure 14. As Fig. [9] distance of the first-order epicyclic over- 
density, yc, estimated with the formalism described in Sec. 12.51 
for the inclined orbit at 8.5 kpc. The corresponding values of Ft, 

V and At can be found in Table [2] The values for the velocity 

V and the acceleration A on the cluster were extracted from our 
computations in time steps of 10 Myr. For comparison the value 
of the reference cluster is also shown. Like on eccentric orbits the 
distance yc starts oscillating about a mean value as a result of 
the acceleration and deceleration due to the galactic disk. The 
mass of the cluster in this plot is fixed to the initial mass of 
20000 Mq. If mass loss was taken into account the mean value 
as well as the amplitude of the oscillations would decrease with 
time. The predictions made within this figure can be compared 
with the A^-body results in Fig. 1131 



This asymmetry may be interpreted as another overdensity 
type, which may be very helpful for determining the orbit 
of an observed cluster like Pal 5. 

The most important point about the studied pericen- 
tre passages is that all observed overdensities are due to 
epicyclic motion of escaping stars or an asymmetric accel- 
eration of the tidal tails. There is no detectable overden- 
sity being obviously the result of a pericentre passage (if we 
neglect the loss of the primordial escapers during the first 
revolution). Hence, even for the most eccentric orbits stud- 
ied here, the conditions of slow escape are well fulfilled for 
a majority of escapers such that epicyclic overdensities can 
be observed. This further supports the assumptions made 
in Sec. 12.51 which will be further investigated in Kiipper, 
Kroupa, Baumgardt & Heggie (in preparation). 



3.5 Milky Way potential, inclined orbit 

We study disk shocks using a cluster on a circular orbit with 
90 degree inclination to the disk. The tidal field which this 
cluster experiences is essentially the same as the one of the 
reference cluster, except for two disk shocks per revolution 
which put additional energy into the cluster according to 
equations [3] & [4] 

In Fig. |3] we can see that the effect of disk 
shocks at 8.5 kpc on the mass evolution of a cluster 
of the given concentration is not significant, in agree- 



ment with results from Vesperini fc Heggie] l|l997f ) and 
iGnedin. Lee fc Ostrikerl (|l999T l. 

Similar as in the reference cluster, the factor d 2 $/dR 2 
(which is important for the calculation of the tidal radius but 
has nothing to do with the shock itself as it only reflects the 
potential change in the radial direction) is as good as con- 
stant within one orbit and just the angular velocity changes 
slightly when the cluster is accelerated or decelerated by the 
disk. This results in the orbit being slightly non-circular due 
to the non-spherical potential. The effective eccentricity is 
only about 0.04, though (see Table [1]). This small non-zero 
eccentricity has a negligible effect on the evolution of the 
cluster. 

The internal evolution of the cluster is also quite similar 
to the reference cluster (Fig. [SJ). Hence, epicyclic overdensi- 
ties are observable in the number distribution of stars along 
the tidal tails, as can be seen in Fig. 1131 where the evolution 
of the tidal tails with time is shown. As the disk periodically 
accelerates and decelerates the cluster quite similarly to an 
eccentric orbit, the overdensities and the overall shape of the 
tails also get slightly compressed and stretched. 

With the formalism which we introduced in Sec. l2.5l and 
the values for the mean galactocentric radius and the mean 
orbital velocity we can predict the behaviour of yc for this 
model. The result is shown in Fig. [14] and is in satisfactory 
agreement with the A^-body results as shown in Fig. 1131 

Important about this model is that no overdensity due 
to a disk shock can be observed within the tails, i.e. there is 
no group of escapers moving along the tidal tail. All over- 
densities are due to epicyclic motion of escaping stars and 
mass loss due to a shock does not happen instantaneously 
with the shock, but happens delayed along the whole orbit 
as described in Section 12.1.31 Furthermore, disk shocks do 
not significantly increase the scatter in escape conditions for 
the epicyclic overdensities to vanish. 

3.6 Orbits at 4.25, 12.75 and 17 kpc 

To see whether the behaviour of the cluster, and especially 
of the tails, changes when the strength and the duration 
of the pericentre passages and disk shocks are varied, we 
investigate all of the above models at three further galac- 
tocentric distances of 4.25, 12.75 and 17 kpc. The initial 
conditions of the test cluster are adjusted to the given tidal 
field, i.e. a Plummer model with a half-mass radius s uch that 
R h/Rtide = 0-2, since Tanikawa fc Fukushigel (|2005l ) as well 
as lGieles fc Baumgardt (2008) have shown that the ratio of 
half-mass and tidal radius significantly determines the evo- 
lution of a cluster. Hence, our setup is a reasonable choice 
for a good comparability. On the other hand, due to the 
fixed ratio of half-mass radius and tidal radius we simul- 
taneously change the two-body relaxation time-scale with 
these adjustments. 

Most globular clusters of t he Milky Way lie below the 
ratios studied here (|Harrislll996h . i.e. are more compact and, 
hence can resist tidal perturbations more easily and have a 
smaller two-body relaxation time. 



3.6.1 Internal evolution 

In Fig. [15] the mass evolution of all clusters at 4.25, 12.75 
and 17 kpc is shown. As mentioned above, since all clusters 
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Figure 13. As Fig. [7] time evolution of the number distribution of stars along the tidal tails in bins of 25 pc for the cluster with 
galactocentric distance of 8.5 kpc on an orbit with 90 deg inclination to the galactic disk (discussed in Sec. 13. 5^ , The first- and second- 
order overdensities can be seen, periodically changing in distance, as the orbit is slightly eccentric. But in mean the evolution is following 
the trend of the reference cluster (cf. Fig.[7Jl. Furthermore, there is no indication that a disk shock creates an ovcrdcnsity by a temporary 
increase in the escape rate. 
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Figure 15. Mass evolution of the clusters at an initial galactocentric distance of 4.25 kpc (left), 12.75 kpc (middle) and 17 kpc (right). 
For comparison the reference cluster is also shown. The increase of importance of tidal influences with decreasing galactocentric radius 
is obvious, where disk shocks gain more influence compared to pericentre passages. 



have the same initial concentration the two-body relaxation 
time decreases for the clusters closer to the galactic centre. 
Therefore, the dissolution of the models on circular orbits is 
stronger the closer the model is to the galactic centre. 

Moreover, an increasing influence of pericentre passages 
for decreasing galactocentric radii can be seen when com- 
paring the models on eccentric orbits at a given apogalactic 
distance with the respective circular model. Models with 
a smaller apogalactic distance go deeper into the bulge at 
perigalacticon and hence are more strongly influenced by the 
tidal-field variations as the angular velocity and hence the 
stress on the cluster grows with decreasing galactocentric 
radius. 

In Fig. [16] the mean stellar mass within the clusters is 



plotted versus the bound mass left in the cluster. An earlier 
increase in mean mass shows that the two-body relaxation 
time scale is short compared to the time scale of dissolu- 
tion of the cluster. Therefore, the probability of observing 
epicyclic overdensities grows for models with an earlier in- 
crease in mean mass. As we can see from the figure, all 
models show an increase in mean stellar mass. Especially at 
4.25 kpc, where the two-body relaxation time is shortest, 
most models show a steep increase in mean stellar mass and 
epicyclic overdensities are therefore expected. 

From our sample just the model at 4.25 kpc with ec- 
centricity of 0.75 dissolves so quickly (Fig. [15} and shows so 
little increase in the mean mass (Fig. I16[) that we would not 
expect observing epicyclic overdensities in this case. 
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Figure 16. The average mass of stars, M/N, versus bound mass, M, for the clusters at an initial galactocentric distance of 4.25 kpc (left), 
12.75 kpc (middle) and 17 kpc (right). The preferential loss of low mass stars is very pronounced in most clusters. Since preferential loss 
is a phenomenon of two-body relaxation, the earlier the increase in M/N with decreasing M the more important is two-body relaxation 
in the cluster. 
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Figure 17. As Fig. [7] time evolution of the number distribution of stars along the tidal tails in bins of 25 pc for the cluster with 
apogalactic distance of 4.25 kpc and an eccentricity of 0.25 (discussed in Sec. I3.6| l. Three to four orders of overdensities can be seen in 
this diagram, periodically changing in distance due to the eccentricity. The clusters at this apogalactic distance have the smallest radii 
and therefore the shortest relaxation times. The periodic influence of the tidal variation does not have a large effect on the cluster, such 
that the escape conditions for the formation of epicyclic overdensities are still well fulfilled for a majority of escapers. 



3.6.2 Evolution of the tidal tails 

As shown in KMH, in the case of a cluster in a constant 
tidal field, i.e. on a circular orbit, epicyclic overdensities are 
generated by stars which escape from the cluster as a result 
of two-body relaxation. The tidal variations studied in this 
work enhance the mass loss of the investigated clusters and 
also the scatter in escape conditions. For most orbital pa- 
rameters two-body relaxation seems to be actively involved 
in the dissolution of the cluster such that the influence of 
the tidal variations can not be too violent. 

We indeed find epicyclic overdensities for all our models. 
In Fig. [17] we show as an example the disk-shock model 



at 4.25 kpc in a time series of the number density along 
the tidal tails. The figure shows two pronounced orders of 
overdensities which oscillate, as the cluster is accelerated 
and decelerated by the galactic disk, and in mean decrease 
in distance as the cluster loses mass. 

Contrary to our expectations, even the most extreme 
cluster in our sample, the one at 4.25 kpc with an eccen- 
tricity of 0.75, is able to grow epicyclic overdensities. Its life 
time is very short; the cluster survives for only about three 
revolutions. Still this is sufficient for most of the escaping 
stars to go through at least one epicyclic loop. 

In Fig. [18] this cluster is shown at t = 140 Myr when it 
has a mass of about 3000 Mq and has completed about 5/4 
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Figure 18. Number distribution of stars along the tidal tails in 
bins of 25 pc (upper panel) and corresponding mean velocities 
of the stars within the bins (with respect to the orbital veloc- 
ity of the cluster, eq. [5} for the cluster with the most extreme 
tidal conditions in our set, i.e. an apogalactic distance of 4.25 
kpc and an eccentricity of 0.75 (discussed in Sec. 13, 6^ . The snap- 
shot was taken at t = 140 Myr when the cluster had a mass of 
about 3000 Mq. Error bars in the upper panel are Poisson errors 
and in the lower panel show the velocity dispersion within the 
corresponding bin (eq. [6}. Notice the different scale on the y-axis 
compared to the other figures, which is necessary as mass loss 
from this cluster is more intense. But even though escape from 
this cluster is happening with a wide range of escape conditions 
epicyclic overdensities are still visible in the tails. 

revolutions about the galaxy (one epicycle of an escaping 
star takes about 1/2 revolution). The mean stellar density 
in the tails in this case is comparatively high as many stars 
are lost from the cluster, such that we had to increase the 
scale on the y-axis. The first-order epicyclic overdensities are 
prominently imprinted. The orbital velocity of stars shows 
no clear variation along the tails except the large-scale vari- 
ation due to the orbit about the galaxy. The velocity dis- 
persion in the corresponding bins are quite large, though. 
This indicates that the scatter in escape conditions is large. 
Nevertheless, there must be a dominant speed and a domi- 
nant radius from which stars escape such that they interfere 
constructively and the observed overdensities can build up. 

In the end, our numerical investigation has convinc- 
ingly shown that we have to anticipate epicyclic overden- 
sities whenever we observe tidal tails. 



4 CONCLUSIONS 

Our numerical investigation has shown that star clusters 
always grow epicyclic overdensities of the kind predicted in 
KMH. We were able to prove this for a large variety of orbits 
and for a test cluster which is comparatively fluffy (large 
ratio of half-mass radius to tidal radius) and hence more 
vulnerable to tidal influences than most star clusters of the 
Milky Way. Only in very extreme conditions, if the cluster 
gets literally disrupted by a single shock, the overdensities 



would not be able to arise, as they need a certain time (about 
one to two orbits) to build up. Since this would imply that 
such a cluster has a very low chance of being observed, we 
can conclude that epicyclic overdensities can be found in any 
star-cluster tidal tail on the sky. 

In KMH we argued that the epicyclic overdensities arise 
if a majority of stars leaves the cluster with evaporative es- 
cape conditions, i.e. with a velocity just slightly above the 
escape velocity when escaping from the cluster through one 
of the two Lagrange points. For clusters on circular orbits 
these conditions are established via two-body relaxation. 
Pericentre passages and disk shocks perturb the evolution 
of a cluster and increase the scatter in escape conditions. 
Nevertheless, we observe epicyclic overdensities for all tested 
clusters. The influence of pericentre passages and disk shocks 
on the escape conditions therefore has to be rather moder- 
ate for the investigated range of orbital parameters. Since 
most clusters of the Milky Way are more compact than our 
test cluster, the average influence of eccentric orbits or disk 
shocks on the escape conditions should be even weaker than 
in our investigation. 

In KMH we showed that the distance from the cluster 
centre to the first epicyclic overdensity, yc, is a multiple of 
the tidal radius. Hence, as the cluster constantly loses mass 
and the tidal radius consequently gets smaller, yc monoton- 
ically decreases with time. Here we showed that this is just 
the case for circular orbits, since for eccentric orbits the pe- 
riodic acceleration and deceleration of the tidal tails leads to 
a stretching and compression of the whole system and there- 
fore periodically increases and decreases yc (e.g. Fig. I10[) . 
At apocentre yc is minimal as the system is maximally com- 
pressed while at pericentre yc reaches its maximum value. 
The strength of the variation in yc depends on the eccen- 
tricity and is a factor of about 7.5 for an eccentricity of 0.75 
(Fig. |9]) . However, on the mean, yc decreases with time as 
the cluster evaporates. 

We furthermore provide an approximation for estimat- 
ing the distance yc in the case of time-variable tidal fields 
(Sec. 12. 5p . It is based on the assumption that the cluster 
experiences a mean tidal field on its orbit, as it is not able 
to react to the tidal variations in time. This assumption also 
suggests that escape from a cluster in a time-dependent tidal 
field happens from the "edge" of the cluster, which is compa- 
rable to the apogalactic tidal radius, rather than from a the- 
oretically determined tidal radius. Kiipper, Kroupa, Baum- 
gardt & Heggie (in preparation) show that the radius of 
the "edge" of a cluster in the given range of orbital parame- 
ters stays approximately constant during one period. Hence, 
the assumptions are justified and the theoretical predictions 
agree well with the iV-body results. 

Moreover, we found further over- and underdensities 
which arise through the non-uniform acceleration of the 
stars within very extended tidal tails; since tails may easily 
span several kpc, different parts may experience very dif- 
ferent forces and have significantly different velocities. This 
leads to a non-uniform stretching and compression of the 
tidal tails causing a large-scale variation of the stellar den- 
sity along the tails (Fig. llip . This effect may be very helpful 
for determining the orbit of an observed cluster like Pal 5. 

Throughout our investigation we were not able to de- 
tect overdensities which could be attributed to a pericentre 
passage or a disk shock, i.e. a bunch of stars which gets 
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unbound during one tidal event and is then moving in a 
group along the tail. Furthermore, we have shown that the 
mass-loss rate is approximately constant for all kinds of or- 
bits when the tidal radius is calculated using eq. 1101 The 
mass-loss rate is larger for more eccentric orbits or orbits 
with disk shocks closer to the galactic centre, but does not 
cause steps in the mass curve, if the pericentre dips are ne- 
glected (Fig. [3}. This is due to a delay in escape of stars 
which get energetically unbound through a tidal event. We 
argue that the stars which are most affect by these events 
orbit the cluster at radii comparable to the apogalactic tidal 
radius of the cluster. The orbital times of those stars within 
the cluster is of the same order as the orbital time of the 
cluster within the galaxy. Escape therefore does not happen 
instantaneously, but the stars leave the cluster delayed with 
respect to the tidal event, where the delay of a specific star 
may last several orbital times (Fig. [8}. As a consequence of 
this, we furthermore suggest to use the apogalactic tidal ra- 
dius for the computation of internal quantities rather than 
the perigalactic tidal radius as has often been done in nu- 
merical studies. 

The Milky-Way globular cluster Palomar 5 is the only 
cluster with prominent, extended tidal tails which, further- 
more, show significant substructure. The most detailed nu- 
merical invest i gation of this cluster has been performed by 
iDehnen et al.l (|2004T ) but they were not able to reproduce 
the substructure within their models. Instead of adding sub- 
structure to the Milky- Way potential like spiral arms or gi- 
ant molecular clouds to increase the time variations of the 
tidal field, or adding dark-matter sub-haloes to the Milky- 
Way potential to perturb the dynamically cold tidal stream, 
as suggested by Dehnen et al., we suggest that most over- 
densities within the tails of Pal 5 are of epicyclic origin or 
are due to an asymmetric acceleration along the tails. 

Another interesting object for t his kind of analysis is 
the recently discovered GD- 1 stre am (|Grillmair fc Dionatosl 
l200fj ; IKoposov. Rix fc Hogdl2009l ). which is thin and thus 
kinematically cold. Even though the progenitor of this 
stream has not been found yet it is most probably ema- 
nating from a globular cluster which may have dissolved by 
now. This stream also shows substructure which may help 
to locate the position of the cluster within the stream since 
the lowest orders of epicyclic overdensities (the overdensities 
closest to the cluster) are always most prominent. 

There have also been stellar overdensities discovered 
without any trace of a tidal stream or an object to which 
they can be attrib uted to, e.g. the Bootes II dwarf spheroidal 
l|Koch et al.ll2009l 'l . With ongoing and future surveys looking 
for stellar overdensities in the Milky- Way halo, the number 
of such objects may well increase. To ease the identification 
of these objects as epicyclic overdensities within tidal tails 
which are themselves below the detection limit, we are going 
to provide a detailed investigation with predictions for ob- 
servers including such factors as radial velocity signatures, 
the stellar mass function, typical sizes and masses, etc., in a 
future contribution. 
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